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NONLINEAR  RESPONSE  OF  BIMODULAR-MATERIAL  PLATES 


J.  N.  Reddy  and  W.  C.  Chao 
Department  of  Engineering  Science  and  Mechanics 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  VA  24061 

Abstract  The  paper  presents  finite  element  results  for 
geometrically  nonlinear  response  of  fiber-reinforced,  single-layer  and 
two-layer  cross-ply  rectangular  plates  constructed  of  materials  which 
have  linear  elastic  properties  in  tension  and  compression  that  are 
different.  A  shear  deformation  theory  of  layered  composite  plates, 
accounting  for  large  rotations  (in  the  von  Karman  sense)  and  the 
bimodular  action,  is  employed  to  analyze  rectangular  plates  made  of  two 
cord-rubber  bimodular  materials.  Numerical  results  for  transverse 
deflection  are  presented  for  simply  supported  plates  under  sinusoidally 
distributed  and  uniformly  distributed  transverse  loads. 

1. _ INTRODUCTION 

The  present  paper  is  a  continuation  of  the  research  by  the  authors 

A 

and  their  colleagues  [1-5]  In  the  analysis  of  bimodular  composite 
plates.  The  previous  investigations  by  the  authors  and  others  were 
based  on  geometrically  linear  theory  of  plates.  The  only  exception  to 
this  statement  Is  provided  by  the  works  of  Kamiya  [6,7],  which  are 
concerned  with  a  clamped  circular  plate,  and  a  simply  supported 
rectangular  plate  under  sinusoidally  distributed  load,  respectively. 

The  present  paper  employs  the  finite  element  developed  In  [2,8]  and  the 
fiber-governed  constitutive  model  of  Bert  [9]  to  Investigate  the 
geometrically  nonlinear  response  of  bimodul ar-materlal  plates.  The 
following  brief  review  of  literature  provides  a  background  for  the 
present  paper. 


Analysis  of  plates  made  of  bimodular  materials  began  with  the  work 
of  Ambartsumyan  [10]  in  1965  (although  Timoshenko  [11]  considered  the 
flexural  stresses  in  such  materials  as  early  as  1941,  Ambartsumyan's 
work  is  credited  for  the  renewed  interest  in  the  analysis  of  bimodular 
materials).  Shapiro  [12]  considered  the  simple  problem  of  a  circular 
plate  under  a  pure  bending  moment  at  its  edge.  Kamiya  [6,7]  analyzed 
the  large-deflection  behavior  of  clamped  circular  plates  using  a  finite 
difference  technique,  and  rectangular  plates  under  sinusoidally 
distributed  load  using  the  Galerkin  method.  In  these  investigations, 
only  bimodular  isotropic  materials  were  considered,  and  the  transverse 
shear  strains  were  omitted.  The  effect  of  thickness  shear  deformation 
was  included  in  the  simple  case  of  cylindrical  bending  by  Kamiya  [13]. 
The  first  analysis  of  bimodular,  anisotropic  materials  is  apparently  due 
to  Jones  and  Morgan  [14],  who  treated  cylindrical  bending  of  a  thin, 
cross-ply  laminate.  In  the  last  couple  of  years,  a  number  of  papers 
dealing  with  the  static  bending  and  free  vibration  of  single-layer  and 
two-layer  cross-ply  plates  have  appeared  [1-5,  15-17].  Most  of  these 
works  are  a  result  of  the  support  of  the  research  by  C.  W.  Bert  at  the 
University  of  Oklahoma  and  the  senior  author  by  the  Office  of  Naval 
Research.  The  significant  contributions  of  this  research  over  previous 
works  are: 

(i)  the  material  of  each  layer  is  both  elastically  and 
thermoelastical ly  orthotropic  and  bimodular; 

(1i)  both  single-layer  orthotropic  and  two-layer  cross-ply 
laminated  plate  and  shell  constructions  were  considered 
using  a  fibre-governed  constitutive  model; 

(ill)  transverse  shear  strains  are  included; 
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(iv)  simply-supported  and  clamped  boundary  conditions  are 

considered,  and  sinusoidal  distribution  as  well  as  uniform 
distribution  of  transverse  load  and  temperature  changes  are 
considered; 

(v)  static,  transient,  and  free  vibration  responses  are  studied; 
and 

(vl)  both  exact  (for  certain  edge  conditions  and  loadings)  and 
finite-element  analyses  are  presented. 

The  present  paper  Investigates  the  large-deflection  (In  the  von 
Karman  sense)  behavior  of  single-layer  orthotropic  and  two-layer  cross- 
ply  plates. 


2.  GOVERNING  EQUATIONS 

Consider  a  plate  constructed  of  a  finite  number  of  uniform¬ 
thickness,  orthotropic,  bimodular-materlal  layers  oriented  arbitrarily 
with  respect  to  the  plate  axes.  The  plate  coordinates  are  taken  such 
that  the  xy-plane  coincides  with  the  midplane  of  the  plate.  Under  the 
assumptions  that  the  layers  remain  linearly  elastic  during  the 
deformation  and  the  generalized  Hooke's  law  is  valid,  and  that  no 
debonding  occurs  between  layers,  one  can  employ  the  equations  governing 
the  shear  deformable  theory  of  layered  composite  plates  [8,18].  Since 
these  equations  are  amply  documented  in  the  works  cited  earlier  (see, 
for  example,  [1-5]),  only  the  strain -displacement  relations  and  the 
equations  of  motion  will  be  repeated  here  to  Indicate  the  nonlinear 
terms  resulting  from  the  von  Karman  theory. 

Assuming  that  the  conditions  of  the  von  Karman  plate  theory  are 
valid,  and  accounting  for  the  transverse  shear  strains,  the  strain- 
displacement  relations  can  be  expressed  In  the  form, 
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au  1  ,a w.2  a*x  _  o 

ei  n  z  W  z  sir  =  e i  z*] 
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aw 
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aw 
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Here  u,  v,  w  are  the  midplane  displacements  along  x,  y,  z  directions; 

and  <|>  and  <|>  are  the  slopes  in  the  xz  and  yz  planes  due  to  bending 
x  y 

only.  In  writing  the  strain-displacement  equations,  it  is  assumed  that 

the  products  of  <|>  ,  «|»  ,  au/ax  and  av/ay  are  neglected.  Since  the 
x  y 

constitutive  relations  are  based  on  the  plane-stress  assumption,  strain 
e3  does  not  come  into  the  equations. 

Neglecting  the  body  moments  and  surface  shearing  forces,  the 
equations  of  equilibrium  (In  the  absence  of  surface  and  body  forces)  can 
be  written  as. 


where  ,  and 


V  Vy-° 

Vx  *  Nx,y  ■  0 

Q.  +  Q-  +  N(N. ,w)  =  q 
1.x  2,y  '  1  ’  ^o 

M.  +  M.  -  0,  *  0 

l,x  6,y  yl 

-  Ox,  »  0 

6,x  2,y  2 

are  the  stress  and  moment  resultants  defined 


(2.2) 


by 


h/2  h/2 

,M1 )  =  /  (l,z)0l  dz  ,  (Oi.Qo)  *  /  (°5 »°4)dz ,  (2.3) 

-h/2  1  2  -h/2  5  * 


and  N(»)  is  the  nonlinear  operator. 
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JV/(w»Ni )  =  5x  (N1  |?)  +  ly  (n6  I?)  +  ax'  ^n6  +  6y  (N2  (2*4) 

Here  a.  (i  =  1,2,4, 5, 6)  denotes  the  stress  components  (a  =  a  , 

1  IX 

°2  '  °y*  °4  =  V  a5  ‘  "xr  ind  °6  ‘  V- 

3.  FINITE-ELEMENT  FORMULATION 

The  finite-element  model  used  in  the  present  study  is  the  same  as 
that  employed  in  [2]  except  for  the  inclusion  of  the  nonlinear  terms. 
The  formulation  Is  not  repeated  here  but  the  steps  involved  in  the 
nonlinear  analysis  of  the  bimodul ar-material  plates  are  pointed  out. 

The  finite-element  model,  for  a  typical  element,  in  the  present  case  is 


of  the  following  form 


mi*}  »  m. 


(3.1) 


where  (a)  denotes  the  column  of  the  nodal  values  of  the  generalized 
displacements.  The  elements  of  the  element  stiffness  matrix  [K]  are 
given  In  Appendix  I. 

Several  comments  are  In  order  on  the  computational  scheme  used  in 
the  present  study.  First  one  should  note  that  the  stiffness  matrix  [K] 
Is  nonlinear  in  that  it  depends  on  the  displacement  vector.  Therefore, 
an  iteration  technique  must  be  used.  On  the  other  hand,  the  calculation 
of  the  stiffness  coefficients  requires  the  knowledge  of  the  neutral 
surface  locations. 


znx  3  “  ^u,x  +  2  <W,x>  ^4x,x 

2ny  “  "  ^v,y  +  ?  ^w,y^  ^y,y’  ^3’2^ 

which  In  turn  depend  on  the  solution  (u,v,w,4»  ,4.  ).  Thus  another 

x  y 

iterative  scheme  is  required  for  the  determination  of  the  neutral 


6 


surface  locations.  The  latter  iterative  scheme  begins  with  assumed 
values  of  znx  and  zny  (say,  znx  =  zny  =  0)  and  then  A^j,  B^j  and  D-jj  are 
computed  using  these  values.  In  general,  the  neutral -surface  locations 
are  not  independent  of  the  position  (x,y),  and  therefore,  the 
expressions  for  plate  stiffnesses,  A-jj,  B-jj  and  D-jj,  also  depend  on  x 
and  y  coordinates.  Since  the  element  stiffness  coefficients  Kjj  are 
evaluated  at  the  Gauss  points,  the  plate  stiffnesses  Afj,  B-jj  and  D-jj 
are  also  evaluated  at  the  Gauss  points  by  using  the  neutral -surface 
positions  computed  at  the  Gauss  points.  After  obtaining  the  generalized 
displacements,  the  neutral -surface  locations  are  recomputed.  Using 
these  new  values  of  znx  and  Zny,  the  stiffnesses  for  the  next  iteration 
are  computed.  This  procedure  is  repeated  until  the  difference  between 
any  two  consecutive  values  of  znx  (and  zny)  differ  by  a  small 
preselected  value  (say  0.1%).  Once  the  convergence  on  the  neutral 
surface  locations  is  achieved,  iteration  on  the  nonlinear  stiffnesses  is 
carried  until  convergence  on  the  displacements  is  achieved. 

4. _ NUMERICAL  RESULTS  AND  DISCUSSION 

In  the  following,  numerical  results  are  presented  for  rectangular 
plates  made  of  two  bimodular  materials:  aramld  cord-rubber  (AR)  and 
polyester  cord-rubber  (PR),  which  are  used  in  automobile  tires.  The 
material  properties  for  these  two  materials  are  given  in  Table  1.  In 
the  present  study,  a  2  x  2  mesh  of  nine-node  Isoparametric  elements  in 
the  quarter  plate  was  used.  The  shear  correction  coefficients  kf  were 
chosen  to  be  5/6.  All  of  the  computations  were  carried  on  an  IBM  3032 
computer  in  double  precision. 
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Table  1  Material  properties  for  aramid  cord-rubber  and  polyester  cord- 
rubber,  unidirectional,  bimodulus  composite  materials. 


Property  Aramid-Rubber  Polyester-Rubber 

Tensile  Compressive  Tensile  Compressive 


En  (GPa) 

3.58 

0.012 

0.617 

0.0369 

E22  (fiPa) 

0.00909 

0.012 

0.008 

0.01 06 

v12 

0.416 

0.205 

0.475 

0.185 

g12  -  Gi3  (GPa) 

0.0037 

0.0037 

0.00262 

0.00267 

G?3  (GPa) 

0.0029 

0.00499 

0.00233 

0.00475 

A  sunmary  of  the  linear  analysis  is  presented  in  Table  2.  The 
effect  of  the  aspect  ratio  (b/a)  and  thickness-to-side  ratio  ( h/a )  on 
nondimensional i zed  center  deflection  (w)  is  apparent  from  the  results  in 
Table  2  (also  see  Figures  1  and  2).  The  effect  of  the  transverse  shear 
deformation  is  to  increase  the  nondimensional i zed  center  deflection  as 
much  as  30%  for  a  side-to-thickness  ratio  of  a/h  =  25. 

The  results  of  geometrically  nonlinear  bending  of  bimodular  plates 
are  discussed  next.  In  order  to  validate  the  present  element  for  the  no¬ 
nlinear  analysis,  first, single-layer  and  two-layer  cross-ply  rectangular 
(ordinary,  jiat  bimodular)  plates  under  uniformly  distributed  loading  and 
simply  supported  boundary  conditions  were  analyzed,  and  the  results  are 
compared  with  the  analytical  (perturbation)  results  of  Chia  [19,20]  in 
Fig.  3.  The  material  properties  of  the  single-layer  plate  (a/h  =  100, 
a/b  =  1)  are  (see  [19]) 

Ej/E2  =  20.0,  G12/E2  =  0.5,  vJ2  *  0.25.  (4.1) 

The  layer  properties  of  the  two-layer  plate  (a/h  =  100,  a/b  =  1)  (see 

[20]) 
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Table  2.  Comparison  of  closed-form  and  finite  element  solutions  of 
single-layer  (0°)  and  two-layer  (0°/90°)  rectangular  plates 
of  aramid-rubber  material (  1  inear  analysis). 


b/a 

h/a 

single-layer 

(0°) 

Two- 

-1 ayer 

0 

O 

0 

o 

W 

-t 

w 

znx/h 

CFS 

FES 

CFS 

FES 

CFS 

FES 

CFS 

FES 

CPT 

0.4317 

0.7124 

- 

0.4281 

- 

0.6960 

- 

0.01 

0.4317 

0.4316 

0.7134 

0.7137 

0.4281 

0.4280 

0.6969 

0.6973 

0.1 

0.4318 

0.4315 

0.8134 

0.8138 

0.4282 

0.4279 

0.7830 

0.7834 

0.6 

0.2 

0.4319 

0.4316 

1.1125 

1.1128 

0.4284 

0.4281 

1.0342 

1 .0345 

0.4 

0.4322 

0.4319 

2.2862 

2.2862 

0.4290 

0.4287 

1.9802 

1.9800 

0.5 

0.4323 

0.4320 

3.1590 

3.1589 

0.4293 

0.4289 

2.6695 

2.6689 

CPT 

0.4420 

- 

1.8671 

- 

0.4383 

- 

1.7734 

- 

0.01 

0.4420 

0.4420 

1.8689 

1.8698 

0.4383 

0.4383 

1.7751 

1.7760 

0.1 

0.4420 

0.4417 

2.0537 

2.0546 

0.4384 

0.4381 

1.9492 

1.9502 

1.0 

0.2 

0.4421 

0.4418 

2.6058 

2.6069 

0.4385 

0.4381 

2.4635 

2.4647 

0.4 

0.4422 

0.4419 

4.7490 

4.7504 

0.4388 

0.4384 

4.3960 

4.3976 

0.5 

0.4423 

0.4420 

6.3228 

6.3243 

0.4389 

0.4386 

5.7764 

5.7782 

CPT 

0.4454 

- 

3.0096 

- 

0.4434 

- 

2.8916 

- 

.01 

0.4454 

0.4454 

3.01 23 

3.0135 

0.4434 

0.4434 

2.8941 

2.8954 

2.0 

0.1 

0.4454 

0.4452 

3.2774 

3.2783 

0.4434 

0.4431 

3.1478 

3.1490 

0.2 

0.4454 

0.4451 

4.0774 

4.0773 

0.4434 

0.4431 

3.9130 

3.9137 

0.4 

0.4454 

0.4451 

7.2381 

7.2337 

0.4435 

0.4431 

6.9234 

6.9220 

0.5 

0.4455 

0.4451 

9.5801 

9.5724 

0.4435 

0.4431 

9.1397 

9.1365 

f  w  =  (wE^h3)102/(qoa4). 
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Figure  3.  Nonlinear  center  deflection  of  single-  and  two-layer  cross-ply 
square  plates  under  uniform  loading  (e|  =  E^,  1*1,2;  simply 
supported  boundary  conditions) 


Ml 
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E1/E2  =  40.0,  G12/E2  =  0.6,  vg  =  0.2.  (4.2) 

The  present  results  are  in  excellent  agreement  with  those  of  Chia 
[19,20].  Next,  an  isotropic,  bimodular,  simply  supported  square  plate 
under  sinusoidally  distributed  loading  was  analyzed  in  an  effort  to  make 
comparisons  with  the  results  of  Kamiya  [7].  The  present  results  are 
compared  with  those  of  Kamiya  in  Fig.  4.  The  present  results  do  not 
agree  with  those  of  Kamiya  [7]  for  all  ratios  of  EC/E*  (including  1). 
Since  our  results  are  validated  for  Ec/E*  =  1  in  the  previous  example, 
one  must  come  to  the  conclusion  that  Kamiya's  [7]  results  are  in  error. 

The  nondimensional ized  center  deflection  (w/h)  versus  the  load 
parameter  ("q  =  /*•  )  are  presented  In  Table  3  for  single-layer 

(0°)  and  two-layer  (0°/90°)  square  plates  (a/h  =  100)  aramid-rubber  and 
polyester-rubber  materials  under  sinusoidal  loading  and  uniform 
loading.  First  note  that  the  response  of  aramid-rubber  plates  is  more 
nonlinear  than  that  of  polyester-rubber  plates.  Also  note  that  the 
deflection  due  to  uniform  loading  is  about  one  and  one-half  times  that 
due  to  sinusoidal  loading  (also  see  Fig.  5). 


Table  3  Nondimensional ized  transverse  deflection  of  single-layer  (0°) 
and  two-layer  (0°/90°)  square  plates  of  aramid-rubber  and 
polyester-rubber  materials 


(a/h) 

=  100,  w 

=  w/h. 

q  ■  qoE22a 

V). 

q 

Aramid- Rubber 

Polyester-Rubber 

single-layer  (0°) 

SL  UDL 

0°/90° 

SL 

single-layer  (0°) 
SL  UDL 

0°/90° 

SL 

Linear 

0.1869 

0.2959 

0.1777 

0.1081 

0.1683 

0.1776 

10 

0.1776 

0.2720 

0.1660 

0.1069 

0.1654 

0.1707 

20 

0.3350 

0.4963 

0.3100 

0.2110 

0.3233 

0.3254 

30 

0.4743 

0.6821 

0.4353 

0.3115 

0.4718 

0.4639 

40 

0.5977 

0.8417 

0.5452 

0.4081 

0.61 06 

0.5879 

50 

0.7083 

0.9798 

0.6443 

0.5006 

0.7399 

0.6997 

Nondimensional ized  center  deflection 
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Figure  4.  Nondimensional ized  center  deflection  versus  the  ratio 
of  longitudinal  modulus  in  tension  to  the  longitudinal 
modulus  in  compression  for  an  isotropic  bimodular,  square  i 

plate  under  sinusoidal  loading  (simply  supported  case).  J 
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Plot  of  nondimensionalized  center  deflection  versus  the  aspect 
ratio  (a/b)  is  shown  in  Fig.  6  for  single-layer  and  two-layer  (0°/90°) 
aramid-rubber  plates  (a/h  =  100)  under  sinusoidal  loading  (for  the  load 
parameter  value  of  10).  The  effect  of  the  aspect  ratio  is,  relatively, 
more  pronounced  in  the  two-layer  plates  than  in  the  single-layer  plates 
(the  effect  is  to  increase  the  deflection,  w  =  w/h). 

Plots  of  the  nondimensionalized  center  deflection  versus  the  side- 
to-thickness  ratio  and  the  load  parameter  are  shown,  respectively,  in 
Fig.  7  and  8  for  aramid-rubber  square  plates  under  sinusoidal  loading. 
From  the  plots  presented  in  Fig.  7  it  is  clear  that  the  effect  of  shear 
deformation  is  more  pronounced  with  increasing  values  of  the  load 
parameter.  This  can  also  be  seen  from  the  load-deflection  curves 
presented  in  Fig.  8. 


5.  SUMMARY  AND  CONCLUSIONS 

Results  of  the  finite-element  analysis  of  the  equations  governing 
the  Timoshenko-type  shear  deformable  theory  that  accounts  for  geometric 
nonlinearities  of  the  von  Karman  plate  theory  are  presented  for  aramid- 
rubber  and  polyester-rubber  bimodular  composite  plates  under  transverse 
loading.  Both  single-layer  and  two-layer  cross-ply  plates  are  analyzed 
under  sinusoidally  and  uniformly  distributed  loads.  The  effect  of  the 
thickness-shear  (which  is  to  increase  the  deflection  w/h  with  increasing 
values  of  slde-to-thickness  ratio,  a/h)  is  more  apparent  for  side-to- 
thickness  ratios  smaller  than  twenty,  and  for  larger  load  parameter 
values.  The  finite-element  analysis  of  ordinary  -  (i.e. ,  not 
bimodular-)  material  plates  whose  elastic  properties  are  taken  to  be  the 
average  of  compressive  and  tensile  properties  listed  in  Table  1  show 


SL=si ngl e-1 ayer ( 0° ) 
TL=two-layer(0/90) 


SL,  a/h=5- 


TL,  a/h=5j 


XSL,  a/h=100 
"TL,  a/h=100 


Figure  7.  Load-deflection  curves  for  square  plates  of  aramid 
rubber  bimodular  material  under  sinusoidal  loading 


linear. 


Figure  8.  Effect  of  side-to-thickness  ratio  on  the  nonlinear 
deflection  (w/h)  of  square  plates  of  aramid-rubber 

bimodular  material  (a/h*100) 


that  the  deflections  predicted  are  about  one-fourth  (for  aramid-rubber 
plates  with  a/h  =  10)  of  those  predicted  using  the  biraodular  properties 
(see  [16]).  Thus  the  effect  of  bimodularity  is  significant  on  the 
response. 

As  pointed  out  in  the  introduction  of  this  paper,  there  is  only  a 
little  to  be  done  in  the  way  of  finite -element  analyses  of  bimodular 
plates.  The  nonlinear  transient  response  of  bimodular  plates  seems  to 
be  the  final  step  in  the  series  of  investigations  based  on  the  fiber- 
governed  constitutive  model.  As  far  as  the  constitutive  models  are 
concerned,  there  is  still  a  need  for  improved  and/or  realistic  models. 

If  a  single  functional  relationship  between  stresses  and  strains  were 
available  for  both  compressive  and  tensile  regions,  the  analysis  would 
be  much  cleaner  (free  of  any  assumption  concerning  the  state  of  stress 
or  strain  in  the  material). 
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APPENDIX  -  I:  ELEMENTS  OF  STIFFNESS  MATRIX 


[K11] 

[K12] 

[K13] 

CK14] 

CK15 

[K21] 

CK22] 

[K23] 

[K24] 

DC25 

CK31] 

[K32] 

[K33] 

CK34] 

[K35 

CK41] 

[K42] 

[K43] 

[K44] 

[K45 

CK51] 

[K52] 

[K53] 

[K54] 

[K55 

The  matrix  coefficients  are  given  by 

*  vl 

[K11]  «  An[SXX]  +  A66[SW], 

[K12]  -  A12[Sxy]  +  Ag6[Sxy]T  =  [K21]T, 

[K13]  =  An[RjX]  +  A12[Rjy]  +  A6fi([Rjy]T  +  [Ryy])  =  \  [K31f , 
CK14]  =  BnCSXX]  +  =  l*41]1, 

[K15]  =  B12[Sxy]  +  Bfi6[SXy]T  =  CK51]T 
[K22]  *  A^l  ♦  A66[SXX], 

[K23]  *  A12[Rjy]T  ♦  A22[Ryy]  +  Ag6([Rxy]  ♦  [R™])  =  \  [K32]T 
[K24]  *  B12[Sxy]T  +  Bg6[Sxy]  *  [K42]T. 

DC25]  -  B2?[Syy]  ♦  Bg6[SXX1  -  [K52]T 


A 


b 

H 

H 

L 


»P  =  A55[s  ]  +  A44[S^ 


[k o3]  =  i  / 


d(>.  3$.  3*.  8<h 

Cl  sx1  sir  *  "e  W-  5T>  *  "2  BT  srl  dx<l 


[K^4]  =  A55[SX0]  =  [kJ3]T, 

[K34]  =  Bn[RjX]  +  B^CR^]1  +  B66([^]  +  [Rjy])  =  2[K43] 

P<f]  =  A44[Sy0]  =  [Kf  ]T, 

[K35]  =  B12[RXy]  +  B^R™]  +  B6g([Rxy]T  ♦  [R™])  =  2Ck|3] 
[K44]  *  Dn[SXX]  ♦  D66CSyyl  +  A55[S], 

[K45]  =  D12[SXy]  +  D6g[SXy]T  =  [K54]T, 

[K55]  =  066CSXX]  ♦  D22[Syy]  +  A44[S]. 


where 


Sir^3TbirdX<ly  ’  t,T]  =  0,X,y  ’  Sij  5  sij’ 

Re 

p  ,  ,.w 

RC  s  /  2  (~6c^  a?,"  dxdy  ’  C.C.n  =  0,x,y, 

Re 


•*1  =  A11  ^l?)2  +  a12 
^2  =  A12  ®2  +  A22 
^6  =  A16  ^a-?)2  +  A26 


,aw,2  .  aw  aw 

(^y)  +  za16  ^  |jy. 

(**,2  ?A  aw- 

'a y>  ^"26  ax  ay’ 


/aw»2  .  aw  aw 
{tiy'  +  ^"66  ax  ay* 


Here  A^j,  Bij,  and  Oij  are  the  plate  stiffnesses. 
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(Aij,Bij  »°1j)  *  ^_h/2  (1,z,z  *Q1jkJidz  (i,J  "  1,2‘6^ 


A1j  =  /  _h/2  Q1jkXdz  OJ  =  4’5) 


(A.4) 


where  h  is  the  total  thickness  of  the  plate,  denotes  the  plane- 

stress  reduced  stiffness  (1,j  refer  to  the  position  In  the  compliance 
matrix;  k  refers  to  the  sign  of  the  fiber-direction  strain:  k  *  1, 
tensile  and  k  =  2,  compressive;  and  i  refers  to  the  layer  numnber) ,  ej 
and  are  the  strains  and  curvatures  associated  with  the  displacements 
in  (2.1),  and  k^  are  the  shear  correction  coefficients  (see  [2-5]). 
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